Observation of a controllable 7r-junction in a 3-terminal Josephson device 
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Recently Baselmans et al. [Nature, 397, 43, (1999)] showed that the direction of the supercur- 
rent in a superconductor/normal/superconductor Josephson junction can be reversed by applying, 
perpendicularly to the supercurrent, a sufficiently large control current between two normal reser- 
voirs. The novel behavior of their 4-terminal device (called a controllable 7r-junction) arises from 
the nonequilibrium electron energy distribution established in the normal wire between the two 
superconductors. We have observed a similar supercurrent reversal in a 3-terminal device, where 
the control current passes from a single normal reservoir into the two superconductors. We show 
theoretically that this behavior, although intuitively less obvious, arises from the same nonequilib- 
rium physics present in the 4-terminal device. Moreover, we argue that the amplitude of the 7r-state 
critical current should be at least as large in the 3-terminal device as in a comparable 4-terminal 
device. 
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When a normal metal is put in contact with one or 
more superconductors, the properties of both materials 
are modified near the interface. The physical phenom- 
ena associated with superconductor (S)/normal (N) sys- 
tems, namely the proximity and Josephson effects, were 
intensely studied in the 1960's and 70's.El Interest in S/N 
systems was rekindled in the 1990's due to the ability 
to fabricate complex structures with submicrometer di- 
mensions. A new, deeper understanding of the pEosimity 
effect on mesoscopic length scales has emerged,Bu con- 
centrating on equilibrium and linear-response physics. 
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FIG. 1. Scanning electron microscope picture of the sam- 
ple, with schematic drawing of the measurement circuit. The 
sample consists of a T-shaped Ag wire with lateral dimen- 
sions of 50 nm x 70 nm, connected to two 70 nm thick Al 
electrodes and one 230 nm thick Ag electrode. 

Nonequilibrium phenomena in S/N systems are now 
taking the spotlight .013 A major discovery was made by 
Baselmans et al.,u who measured a 4-terminal diffusive 
metal S/N/S Josephson device with a cross shape. Two 



opposing ends of the cross were connected to S electrodes, 
while the other two were connected to N reservoirs be- 
tween which a control current was passed. Baselmans 
et al. found that, at high control current, in samples 
with the normal reservoirs sufficiently close together, the 
sign of the Josephson supercurrent between the S elec- 
trodes reversed direction. The current-phase relationship 
under such conditions becomes I s (4>) = I c sin(r/> + tt), 
where I c is the (positive) critical supercurrent, rather 
than the usual Josephson relationship I s {4>) — I c sin(0), 
hence the device is called a 7r-junction. Such-a device has 
been used to make a controllable 7r-SQUID.EI The expla- 
nation, of the nonequilibrium 7r-junction consists of two 
parts. a First, the supercurrent can be decomposed into 
an energy-dependent "spectral supercurrent" je, which 
is an equilibrium property determined by the sample ge- 
ometry and resistance as well as the phase difference <f> 
between the two S electrodes, je is an odd function of 
energy, and exhibits damped oscillations on an energy 
scale comparable to the Thouless energy of the sample, 
E t h = hD/L 2 , with D the diffusion constant in the wire 
and L the length between the superconductors. Sec- 
ond, the total supercurrent is determined by the occu- 
pation of the supercurrent-carrying states, given by the 
antisymmetric part of the quasiparticle distribution func- 
tion f(E) in the normal region of the junction describ- 
ing the pairs of quasiparticles (E > Ep) and quasiholes 
[E < Ep). Under nonequilibrium conditions, f(E) can 
be made to have a staircase shape, with-stcps appearing 
at the voltages of the normal reservoirs JIj The staircase 
shape of f(E) excludes the low-energy contribution of 
je from the supercurrent. When the control voltage ap- 
proaches the energy where je changes sign, the supercur- 
rent changes its sign relative to the equilibrium situation. 
In contrast to the 7r-junction behavior, smearing of the 
distribution function by electron heating or raising the 
sample temperature simply causes the supercurrent to 
decrease toward zero without ever changing sign. 

The sample shown in Fig. 1 consists of a T-shaped Ag 
wire, 70 nm wide and 50 nm thick, connected to two S 
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electrodes (70 nm of Al) and one N reservoir (230 nm of 
Ag). The distance between S electrodes is 1.1 /im, while 
the distance from the top of the "T" to the N reservoir 
is 4.5 /im. The phase coherence length in similarly 
prepared Ag wires is several micrometers at sub-Kelvin 
temperatures, hence we expect to observe a substantial 
Josephson effect between the two S electrodes. The sam- 
ple was fabricated using one electron-beam and two op- 
tical lithography steps. The T-shaped Ag wire was fabri- 
cated first, followed by the thick Ag reservoir, and finally 
the Al electrodes. A gentle ion mill of the exposed ends of 
the Ag wire preceded the evaporation of the Al electrodes 
to enhance the transparency of the Ag/Al interfaces. The 
sample was immersed in the mixing chamber of a dilution 
refrigerator with filtered electrical leads. 

The transport properties of the sample were deter- 
mined initially by measuring the V vs. I characteris- 
tics between pairs of electrodes. The V-I curve between 
S electrodes shows the standard Josephson junction be- 
havior with a critical current of 0.7 /iA at 38 mK. The 
V-I curve between the N electrode and either S elec- 
trode exhibits a change in slope at a current approxi- 
mately equal to twice the critical current. This behavior 
is due to the superposition of opposite-flowing quasipar- 
ticle current and supercurrent in the dangling arm, as ob- 
served recently by Shaikhaidarov et al.tLil. For the sample 
shown in Fig. 1, the left and right arms have resistances of 
Ri = 7.0 ft and i?2 = 9.1 f2, respectively, while the base 
of the T has a resistance of Rq = 36 17. From these val- 
ues and the sample geometry, we deduce that about half 
of the 16.1 fl S-S resistance comes from the uncovered 
part of the Ag wire, and the other half from the Al/Ag 
interfaces and part of the Ag wire extending under the 
Al electrodes. 




I (MA) 

FIG. 2. A subset of Vsns vs. / curves measured across the 
S/N/S Josephson junction, for different values of the current 
injected from the normal reservoir. From bottom to top, the 
injected currents Iinj are in (j,A: 0.53, 0.70, 1.01, 1.23, 1.89, 
2.18, 3.15. The curves are offset for clarity. 



The measurement circuit for the nonequilibrium injec- 
tion experiment is shown schematically in Fig. 1. A dc 
current I m j is injected from the normal electrode to one 
of the superconducting electrodes. Simultaneously, the 
V-I curve between the two superconducting electrodes is 
measured in a 4-probe configuration. Figure 2 shows a 
subset of V-I curves for different values of ii n j , and is the 
central result of this paper. The critical current of the 
S/N/S junction decreases rapidly with increasing injec- 
tion current. When I m ^ = 1.0 /iA, the critical current is 
below our measurement threshold. Upon further increase 
of Ijnj, the critical current increases again, and finally 
disappears when I ln] > 3 /iA. In Figure 3, we plot I c vs. 
Vn at three different temperatures, where Vn — iijv-^inj 
is the voltage of the normal reservoir with respect to the 
superconductors, and Rn = Ro + (Ri 1 + R^ 1 ) -1 = 40 f2. 
In the figure we intentionally plot I c < after it falls to 
zero, to emphasize that the junction has entered the 'V 
state.Ej Our interpretation of the data is consistent with 
the assumption that, for fixed <f>, I s is a smooth func- 
tion of Vn with a continuous first derivative. It is also 
consistent with the experiment of Baselmans et al.p who 
confirmed the existence of the 'V state by measuring 
the resistance of the normal wire as a function of the su- 
percurrent, hence the phase difference 0, between the S 
electrodes. At zero supercurrent, their wire resistance ex- 
hibits a local minimum in the usual "0" state and a local 
maximum in the 'V state due to the proximity effect. 
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FIG. 3. Critical current of the Josephson junction vs. volt- 
age of the normal reservoir at T=38 (□), 96 (A) and 200 mK 
(o). I c is shown as negative for Vn > 40/iV to symbolize the 
appearance of the 7r-junction. Inset: Critical current vs. tem- 
perature at Vn = 0. The lines are the theoretical calculations 
discussed in the text. 

The significant difference between our experiment and 
that of Baselmans et al., aside from the reduction from 
4 terminals to 3, is the presence in our sample of a dissi- 
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pative quasiparticle current in the sample arms that si- 
multaneously carry the supercurrent. In the Baselmans 
experiment, the control voltages of the two normal reser- 
voirs were set to values ±V/v with respect to the super- 
conductors, so that the electrical potential was zero ev- 
erywhere along the wire connecting the two superconduc- 
tors. To compare our experiment with theirs, we must 
understand the influence of the dissipative current on the 
supercurrent in our sample. We use the quasiclassical for- 
malism in real time, which was originally developed fpc 
nonequilibrium phenomena in massive superconductorstJ 
but also adapted and successfully applied to mesoscopic 
proximity systems, as reviewed e.g. in refs. 3 and 14. 

For the present paper, we are concerned primarily with 
the supercurrent 
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where cat and A are the conductance and cross-section 
of the normal wire, je is the spectral supercurrent dis- 
cussed earlier, and fL{E) = f(—E) — f(E) is the anti- 
symmetric part of the electron energy distribution func- 
tion. With the chemical potential of the superconduc- 
tors taken to be zero, the symmetric distribution function 
fr(E) = 1 — f(E) — f(—E) describes charge imbalance, 
while Jl(E) describes energy or heat in the conduction 
electron system. 

To calculate the supercurrent, first one must solve the 
Usadel equation for the retarded and advanced Green's 
functions. Those contain all information about energy- 
dependent properties of the sample, including the func- 
tion 2e- To find f^(E), one must then solve the Keldysh 
component of the Usadel equation, which takes the form 
of consepation laws for the spectral charge and heat 
currents.li-3 When j B ^ 0, the two kinetic equations are 
coupled, and lead to complicated spatial and energy de- 
pendences of /l(E) and fr(E) in the arms of the sample 
between the superconductors. A major simplification oc- 
curs in the arm of the sample connected to the normal 
reservoir: je = there since the superconducting phase 
is constant along that arm. For voltages and temper 
atures small compared to A the heat current is zero,cJ 
hence Jl(E) is constant along that arm and takes on 
the (equilibrium) value it has in the N reservoir: f L — 
{l/2){tanh[(E + eV N )/2k B T}+ta,nh[(E-eV N )/2k B T]}. 
Since the total charge current is conserved along the 
two sample arms connecting the superconductors, we can 
evaluate it anywhere in those arms. At the central point, 
the dissipative currents diverted into the two arms can- 
cel and we can find the supercurrent from Eq. (1) using 
the expression for f L (E) given above, without integrat- 
ing the kinetic equations. We need only to evaluate je 
at the central point by solving the equilibrium Usadel 
equation for our sample geometry. _ 

As an extension of previous work,0 we have solved the 
retarded Usadel equation taking into account the influ- 
ence of the lead to the normal reservoir and the finite 



interface resistancesB The normal reservoir induces ex- 
tra decoherence into the structure, decreasing the mag- 
nitude of the observed supercurrent. -We find that the 
full gap in the spectral supercurrento becomes a pseu- 
dogap and that the amplitude of the maximum of j'e is 
strongly reduced (although the total supercurrent is re- 
duced by only 20% at 40 mK) . Our fit to the equilibrium 
data of critical supercurrent vs. temperature is shown in 
the inset to Figure 3. To fit the temperature dependence, 
the Thouless energy was adjusted to be £>rh = 3.5 /LteV, 
which corresponds to a distance L — 1.7 /im between 
the superconducting electrodes - larger than the actual 
distance as a result of the silver wire penetration un- 
der the aluminium reservoirs and of the finite contact 
resistances. Surprisingly, the magnitude of the calcu- 
lated cricital current had to be reduced by a factor 1.7 to 
match the experimental data, possibly due tq-the rather 
high S/N interface resistances in this sampleO 

If we now calculate the nonequilibrium data of I c vs. 
Vn using the equilibrium form for fl in the normal reser- 
voir, we find that the calculation overestimates the criti- 
cal current in the 'V state by a large factor, and predicts 
too small a voltage at which the supercurrent changes 
sign. This failure results from neglecting inelastic colli- 
sions inside the wire and electron heating in the normal 
reservoir. Based on our previous measur emen ts of f(E) 
in nonequilibrium mesoscopic metal wires ,l 10 H la l we can es- 
timate the contributions of both inelastic scattering and 
reservoir heating to the rounding of f(E) in our sam- 
ple. Inelastic scattering in similar Ag wires was well de- 
scribed within the framework of the Boltzmann equation 
using an electron-electron interaction kernel in agreement 
with the theoretical form K{E) = K 3 / 2 E~ 3 / 2 , but with a 

prefactor K 3 / 2 ~ 0.5 ns^meV -1 / 2 , about 5 times larger 
than predicted by theory. Heating of the normal reser- 
voir can be estimated using the Wiedemann-Franz law 
and a simplified model of electron-phonon scattering in 
the reservoir .E§E3 The temperatur e of the ele ctrons in 
the reservoir is given by T e ff = \/T 2 + b 2 V^ where b 2 
is proportional to the ratio .of the reservoir sheet resis- 
tance to the wire resistance.£3 From our sample paraiit. 
eters and previous measurements of similar samples ,E2I 
we estimate b w 1 K/mV. Using these values of K 3 / 2 
and b, we have calculated f(E) and thereby I c (Vn) in 
our sample by solving the Boltzmann equation with the 
correct boundary conditions at the S/N interfaces but 
neglecting proximity effect in the bulk of the wire. The 
result of that calculation does not fit the data shown in 
Fig. 3. A much larger value of K 3 / 2 — 3ns _1 meV -1 ^ 2 
provides a reasonable fit, but leaves us without a plausi- 
ble explanation for the enhanced electron-electron inter- 
actions. An alternative approach is to use an interaction 
kernel of the form K(E) = K 2 E~ 2 , which describes sam- 
ples containing dilute magnetic impurities BBcl With the 
value K 2 = 0.55 ns _1 , corresponding to a magnetic impu- 
rity concentration of about 0. 1 ppm, we obtain the solid 
curves shown in Fig. 3, which fit the data well at voltages 
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up to the crossover to the it junction. Adding a reason- 
able K 3 / 2 term to K(E) improves the fit only slightly at 
higher voltages. The magnetic impurity concentration 
of 0.1 ppm is plausible, and will limit to about 5 /im 
near the Kondo temperature - still much larger than the 
distance between the two superconducting electrodes. 

The rather poor fit to the data at high voltages may 
reflect the fact that the magnitude of I c in the 7r state 
depends on a delicate balance between the positive and 
negative parts of je, weighted by the precise shape of 
f(E). Fig. 4 shows f(E) for V N = 50 /iV, near the 
maximum 7r junction I c . By eye f(E) looks nearly like 
a hot Fermi-Dirac function, but the dashed line in the 
figure shows that it is not. If the sample were shorter, 
so that f(E) maintained the staircase structure of the 
dotted line in the figure, the it junction I c would be much 
larger. 
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FIG. 4. Left: Solid line: distribution function f(E) used 
to calculate the Josephson junction current in the n state at 
Vn = 50flV and T = 38 mK. Dotted line: f(E) taking into ac- 
count only reservoir heating but not energy exchange. Dashed 
line: hot Fermi-Dirac distribution. Right: Numerically cal- 
culated je (multiplied by the prefactor on A), at the central 
point of the sample, shown only for E > 

Fig. 4 also reveals the difference between our 3- 
tcrminal experiment and the 4-terminal experiment of 
Baselmans et al. In our sample the electrical potential is 
nonzero at the central point, since the injection current 
flows into both S electrodes. Hence f(E = 0)^^ at the 
central point, unlike in Baselmans' sample. (The devia- 
tion from 1/2 is small, since the vertical arm of our sam- 
ple is much longer than the horizontal arms.) Since the 
available phase space for quasiparticle energy exchange 
decreases as f(E) deviates from 1/2, the 3-terminal ge- 
ometry should be favorable for maximizing I c in the 7r 
state. A direct measurement of this subtle effect could 
be made in a 4-terminal sample. Biasing the two nor- 
mal reservoirs at the same potential Vn, rather than at 
asymmetric voltages ±Vat, would result in a current flow 
pattern and distribution functions essentially equivalent 
to those in our 3-terminal experiment. A comparison 
of the values of I c in the ir state under symmetric bias 
(Vn,Vn) and antisymmetric bias (V/v, — V/v) might re- 
veal a subtle difference in the smearing of f(E). We plan 
to explore this comparison experimentally. 
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